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(54) fluid measurement 

(57) In measuring the amount of fluid 
contained in a tank (10), two values 
L,. Lj related to the amount of fluid 
present are calculated* from 
information provided respectively by a 
level gauge 1 1 and a flowmeter 12. A 
calculating unit 1 4 resets the second 
value when the difference between the 
two values exceeds a given threshold. 
Also, a factor applied to the flowmeter 
output is modified depending on the 
difference existing between the two 
values. The display 1 7 displays the 
second value L-. 
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SPECIFICATION 
Fluid mMsur»m«nt 

The present invention relates to the 
measurement of the amount of fluid contained In 
5 a tank, particulariy the amount of fuel present in 
the tank of a motor vehicle. 

(n the motor car field, the Important 
development has been witnessed for some years 
of apparatus known as "car computers" which are 
1 0 intended to supply the driver with a targe 
number of items of Information, Including 
average instantaneous consumption, cruising 
range, etc. Nevertheless, in order to ensure a 
correct operation of these apparatus and to 
1 5 facilitate the observation of these, it would be 
desirable to have available measuring systems 
which are precise and deliver stable Information. 

The more conventional devices for measuring 
the level of the fuel in the tank, such as devices 
20 with a float for example, cannot give satisfaction. 
Indeed, apart from problems associated with their 
complexity and poor reliability, these devices are 
very sensitive to the fluctuations in fuel level 
owing to disturbances, e.g. due to accelerations of 
25 the vehicle. 

An attempt has been made to overcome this 
disadvantage by inserting the float in a tubular 
element which is disposed verticaiiy and In 
communication with the tank through a calibrated 
30 orifice. Nevertheless, this arrangement proves 
relatively complex and It Is found, in practice, that 
it doss not provide sufficiently stable information. 

Other devices have been proposed but these 
have been set aside particularly because of their 
35 high cost and technical complexity. 

In accordance with this Invention 
there is provided o method of measuring the 
amount of fluid contained in a tank, comprising: 

(a) measuring the level of the fluid in the tenk 
40 by means of a gauge. 

(b) determining a first value L,(n) related to the 
amount of fluid present in the tank from a value 
representing the fluid level provided at a given 
moment (n) by the gauge and from a first value 

45 {L^(n— 1)) previously determined. 

(c) measuring the amount of fluid withdrawn 
from the tank per unit of time, 

(d) determining a second value related to the 
amount of fluid present in the tank from the 

50 measured flow rate, epplying thereto a given 
factor (K); 

(e) comparing the values L,, L^. related to the 
amount of fluid, as determined In steps (b) and 

(d); 

55 (f) resetting the second value to equal the 
first value L, each time the difference between 
these two values L,. exceeds a given threshold 
X. said second value (L,) then representing the 
amount of fluid; 

60 (g) modifying the factor K applied to the 
measured flow rate for the calculation of said 
second value (L,) in the step (d) depending on the 
difference existing between the first vaiue and 
the second value L^: and 
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65 (h) repeating the steps la) to (g). 

TIte embodiment to be described herein is 
simple, robust, reliable and economical, white 
providing preci&e and stable information. In 
parucular it makes use of a gauge and preferably 

70 3 flowmeter, which devices are completely 
conventional in themselves. This embodiment 
enables the contents of the tanic to be known in a 
precise and stable manner to supply the user of 
the car computer with stable Information. 

75 Other characteristics and advantages of the 
present invention will ba apparent on reading the 
following detailed description with reference to 
the accompanying drawings, given by way of non- 
limiting example, and in which: 

80 Figure t is a block diagram of a measuring 
device according to the present invention; 

Figure 2 is a diagrammatic flow chart 
illustrating the operation of the device according 
to the present invention; 

85 Figure 3 represents curves which illustrate the 
various values determined by the device 
according to the present invention: and 

Figures 4j and 4b Illustrate the vaiue of the 
amount of fluid contained in the tank, as 

90 calculated by the device according to the present 
invention, in two particular cases. 

As illustrated diagrammatlcally in Figure 1 . the 
measuring device according to the present 
invention comprises a gauge t ) to measure the 

95 level of the fluid, or liquid, contained in the tank 
10 of a vehicle, as well as a flowmeter 12 which 
measures tlie amount of fluid withdrawn from the 
tank per unit of time. 

The Information delivered by the gauge 1 1 is 
1 0O applied to a calculadng unit 1 4 by means of an 

aiialoyue/digital converter 1 3. Of course, it Is easy 
to determine the number of litres of liquid 
contained in the tank 10, on the basis of the 
information delivered by the gauge 1 1 . if the 
1 05 dimensions and the shape of the tank 1 0 are 
known. 

in a similar manner, the pulses delivered by the 
flowmeter 12 and representative of the flow of 
fuel, are applied to the calculating unit 14. There 
1 1 0 again. It is possible to determine the number of 
litres of fuel contained in the tank 1 0 on the basis 
of the information supplied by the flowmeter, if 
the number of litres originally present in the tank 
1 0 Is known. 

115 The measuring device according to the present 
invention likewise comprises two registers of 
stores 1 5 and 1 6. connected to the calculating 
unit 1 4. T^ie first register 1 5 is adapted to store 
said first valuu (L,) of the amount of fluid present 

1 20 in the tank, determined by means belonging to 
the calculating unit 1 4. depending on the 
Information representative of the fluid level, 
delivered by the gauge 1 1 . weighting this 
information on the basis of a first value previously 

125 determined. 

The second register 1 6 is adapted to store a 
second vaiue (L,) of the amount of fluid present in 
the tank, determined by means likewise belonging 
to the calculating unit 14. depending on the 
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Information (dKrepresentatlve of the flow of fluid, 
delivered by the flowmeter 1 2, by applying a 
given factor (K) to this Information. 

As will be apparent on reading the following 
5 description, the calculating unit 14 likewise 

comprises means for comparing the aforesaid two 
values (L,, L,) of the amount of fluid and means 
for resetting the second value (L.) when the ^ 
difference between these values exceeds a given 

10 threshold, said resetting means making the 
second value (L,) equal to the first (L,) in such a 
manner that the second value represents the 
amount of fluid present In the tank. On the other 
hand, the calculating unit 14 comprises means for 

1 5 modifying the factor (K) applied to the information 
(d) delivered by the flowmeter 1 2 for the 
calculation of the second value, depending on the 
difference existing between the first (L^ and the 
second value (L,). 

20 Finally, as illustrated In Figure 1 , the device 
likewise comprises meens 1 7 for displaying the 
second value representative of the amount of 
fluid present in the tank. 

The method of measuring according to the 

25 present invention will now be described with 
reference to the flow chart Illustrated In Figure 2. 

When voltage is first applied and when the 
device is reset to zero, illustrated 
diagrammattcelly by the first step 50: 

30 "Lo=L,=L,". a value Lg representative of the 
amount of fluid contained in the tank and 
calculated by the calculating unit 14 on the basis 
of the information delivered by the gauge 1 1 Is 
transferred into the registers 1 5 and 1 6. 

35 As illustrated by the step 5 1 . the calculating 
unit 1 4 determines a first value L,{n) of the 
emount of fluid present in the tank depending on 
the information representative of the fluid level 
delivered at a given moment {n) by the gauge, 

40 weighting this information on the basis of a first 
value (n— 1 ) determined previously. More 
precisely, as illustrated in Figure 2. the calculating 
unit 1 4 determines said first value L, on the basis 
of the following formula: 



45 



(p-1) L„ 
L,(n-1)+ 



in which: , ^ . 

L,tn) represents the first value L, at the 

moment n. . . . 

L,(n— 1 ) represents the first value L, at the 

50 moment (n—1), 

L« represents a datum proportional to the 

Information representative of the fluid level 

delivered by the gauge, 
p represents any whole number. 

SS Said first value L,. thus determined, is stor^ In 
the register 1 5 and reactualized sequentially, for 
example every second. 

By way of non-limiting example, p may oe 
equal to 1 6 which gives a response time 



60 equivalent to 1 6 seconds and causes a lag of 0.2 
litre for flows of 50 litres/hour. 

As Illustrated in Figure 3. in which there Is 
represented a first curve representative of the 
amount of fluid contained in the tank, determined 
65 on the basis of the information delivered by the 
gauge, as well as a curve representative of said 
first value, this smoothing effected by the 
calculating unit 1 4 eliminates the rapid variations 
In the values recorded by the gauge 1 1 and due. 
70 In particular, to the displacements and 
accelerations of the vehicle. 

At the step 52. the calculating unit 14 
determines a second value of the amount of 
fluid present in the tank depending on the 
75 information d representative of the flow of fluid, 
delivered by the flowmeter 1 2. applying a given 
factor K to this Information d. 

More precisely, the second value of the 
amount of fluid present in the tank is determined 
80 by the calculating unit 14 on the basis of the 
following formula: 

L.Jn)«Lj(n-1 )-ICd. 

in which; 

Ljin) represents the second value at the 
85 moment n. 

tjtn-l ) represents the second value L, at the 

moment (n— 1 ). 
d represents the information delivered by the 

flowmeter 12. representative of the flow of 

90 fluid. 

K represents said factor. 

At the step 53, the values L, and 
determined at the preceding steps 51 and 52 are 
compared. When the absolute value of the 
95 difference between and is greater than a 
given value X. the second value is reset in step 
54, making it equal to the first value L,. 

Said second value is then representative of 
the amount of fluid contained in the tank 1 0 and 
1 00 this second value is displayed, in- step 55, by 
the means 17. The method according to the 
present Invention Is then resumed before the step 
51 in such a manner as to bring said first L, and 
second values representative of the 
1 05 amount of fluid present in the tank 1 0 
constantly up to date. 

Naturally, if the difference between the two 
values and is less than said constant X Istep 
53), the value of Lj determined directly is 
1 10 displayed by .the moans 1 7. 

By way of non-limiting example, the value X 
may be equal to 



1+- 



20 



Naturally, this value X may be selected differently. 
1 1 5 Furthermore, as previously Indicated, the 
method according to the present invention 
consists in modifying the factor K applied to the 
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information d delivered by the flowmeter 1 2, for 
the calculation of said second value L,, depending 
on the difference existing between the first value 
and the second value L,. 
5 According to a first variant, said modification 
effected on the factor K applied to the Information 
d delivered by the flowmeter 1 2 is affected as 
fellows: 

if the first value is lower than the second 
1 0 value L,. the factor K is taken =aK nominal, a 
relationship in which a is greater than 1 and. for 
example, a^^.3. 

whereas if the first value L, is higher than the 
second value L,. the factor K is taken ^fi K 
\ 5 nominal, a relationship in which fi is less than 1 , 
for example 

Such a amoothing preventa any rise of the 
second value L, and always tenda to bring the 
second value closer to the first L,, on the 
20 average. Of course, the values a and ^ can be 
easily determined for this purpose. 

According to a second variant, the modification 
effected on the factor K applied to the information 
d delivered by the fiowmetar 12 Is carried out by 
25 taking the factor K=sK nominal 



a relationship in which: 

L, represents the second value and 
represents the first value. 
30 According to a third variant, the modification 
effected on the factor K applied to the information 
d delivered by the flowmeter 1 2 Is carried out by 
taking the factor K«sK nominal 

'•2 

L, 

35 a relationship in which: 

Lj represents the second value and 

represents the first value. 
Finally, according to a fourth variant, the 
modification effected on the factor K applied to 
40 the information d delh/ered by the flowmeter 1 2 is 
cam'ed out by taking the factor KssK nominal 

fT+f(Lj--t,)] 

a relationship in which: 

f (L^— L<|) represents a function having the 
45 difference existing between the second value 
and the first value as a variable. 
The contents of the register 1 6 are taken into 
account for the display, by the means 1 7, of the 
contents of the tank 1 0 and for the calculation of 
50 the cruising range. 

In order to avoid the problems associated with 
the stopping of the engine and the opening of the 
electric feed circuit, it is possible to envisage 
transferring the contents of the registers 1 5 end 
55 1 6 to permanent stores. Nevertheless, it is 
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possible to provide that when voltage is again 
applied to the circuit, the first measurement of Lg. 
determined on the basis of the information 
delivered by the gauge 11 , Is transferred to the 

60 registers 1 6 and 1 6. as an initial value of and 
Lj. Such an arrangement renders it possible, in 
particular, to take into account, in a simple and 
rapid mannsr, tl^e suddan variations in the fuel 
level in the tank 1 0 which appear, in particular, 

■65 whan said tank 10 is filled. 

As illustrated in Figure 3. in which the curve 
appears, if the gauge 1 1 is placed at a point 
where the level varies little with the inclination of 
the car, even if the fuel level varies with the 

70 waves, the display effected by the means 1 7 is 
effected with the precision of the gauge 1 1 and in 
a continually decreasing manner which is 
theiefore not disturbing to the driver. 

On the other hand, it appears that the various 

75 breakdowns liable to be encountered with the 
device according to the present invention can be 
simply traced. This if the level gauge 1 1 should be 
cut off, the indication delivered by the display 
rnaans 1 7 would oscillate between two values 

80 corresponding substantially to the maximum level 
of the tank 10. Such a phenomenon is due to the 
fact that at the start, the gauge 1 1 Indicates that 
the tank 1 0 Is full. This Indication Is taken into 
account in the registers 1 5 and 1 6. The contents 

85 of the register 1 6 are brought up to date on the 
basis of the pulses delivered by the flowmeter 12. 
therefore the information delivered by the display 
means 1 7 tend to diminish substantially but is 
reset to the maximum value as soon as the 

90 difference between the registers 1 5 and 16 
exceeds ihe aforesaid value X. In a similar 
manner, if the gauge 1 1 Is short-circuited, the 
display means 1 7 indicate that the tank 10 is 
empty. Such a phenomenon is due to the 

95 registers 1 5 and 1 6 initially taking into account 
the information "tank empty" delivered by the 
gauge 1 1 . 

As illustrated diagrammatically in Figure 4a, 
when the flowmeter 1 2 is cut off or short- 

1 00 circuited, the infonnation delivered by the display 
means 1 7 decreases in steps. Such a 
phenomenon is due to the fact that the contents 
of the register 1 6 are not continuously brought up 
to date on the basis of t^e information delivered 

1 05 by the Flowmeter 1 2. because this is not 

delivering any pulses. The register 1 6 is therefore 
only brought up to date, in steps, when the 
difference existing between the register 1 5 and 
the register 1 6 exceeds the aforesaid value X. 

110 Finally, as apparent in Figure 4b. when the 
flowmeter 12 delivers too great a number of 
pulses per cm^ of luel withdrawn from the tank 
1 0, the information delivered by the display 
means 1 7 appears in the form of saw teeth. Such 

115 a phenomenon is due to the fact that the 

contents of the register 1 6 are reduced much 
more rapidly than the contents of the register 1 5. 
In such a case, when the difference between the 
contents of the register 1 5 and the contents of 

1 20 the register 1 6 exceeds the aforesaid value X. the 
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bringing up to date off the register 1 6 always 
tends to raise the value of ihis. 

The flowmeter 12 may be replaced by any 
ether suitable means capable of delivering 
5 information representative of the amount of fluid 
or fuel withdrawn from the tank per unit of time. 
By way of example, the flowmeter 1 2 could be 
replaced by an element sensitive to the number of 
revolutions of the engine or to the admission 
10 suction. 

Claims 

1 . A method of measuring the amount of fluid 
contained in a tank, comprising: 

(a) measuring the level of the fluid in the tank 
y 5 by means of a gauge. 

(b) determining a first value L,(n) related to the 
amount of fluid present in the tank from a value 
representing the fluid hsvel provided at a given 
moment (n) by the gauge and from a first value 

20 IL,(n— 1 )) previously determined. 

(c) measuring the amount of fluid withdrawn 
from the tank per unit of time, 

(d) determining a second value related to the 
amount of fluid present in the tank from the 

25 measured flow rate, applying thereto a given 
factor (K): 

(e) comparing the values L,. L,. related to the 
amount of fluid, as determined in steps (b) and 

id): 

30 (f) resetting the second value L, to equal the 
first value each time the difference between 
these two values L,. exceeds a given threshold 
X. said second value (L,) then representing the 
amount of fluid; 

35 (g) modifying the factor K applied co the 
measured flow rate for the calculation of said 
second value (L,) in the step (d) depending on the 
difference existing between the first value and 
the second value L,: and 

40 (h) repeating the steps (a) to (g). 

2. A method as claimed in Claim 1 . in which for 
step fc) a flowmeter isused. 

3. A method as claimed in Claim 1 or 2. in 
which said first value L, is determined in step (b) 

45 according to the equation: 

(p-1) ^ 
L,(n)= L,(n-1)+ 

P P 

in which: 

L,(n) represents the first value at the 
moment (n). 

50 Lf in— 1 ) represents the first value L, at the 

moment (n— 1 ), 
Lg represents a datum proportional to the value 

representative of th«.* fiaid level delivered by 

the gauge 111). and 
55 p represents any whole number. 

4. A method as claimed in any one of Claims 1 
to 3. in which the second value related to the 
amount of fluid present in the tank is determined 
according to the equation: 



60 



L,(ni=L,(n-l)-Kd. 



in which: 

Ljin) represents the second value L2 at the 

moment In). 
L,(n* 1 ) represents the second value ILJi at the 
55 moment (n— 1 ). and 

(/represents the measured flow rate. 

5. A method as claimed in Claim 4. in which 
the modification effected in step Ig) of the factor 
K is carried out as follows: 

70 if the first value L, is lower than the second 
value L, the factor K is taken as a K nominal, a 
relationship in which a is greater than 1 , 

whereas if the first value is higher than the 
second value L, the factor K is taken as /i K 

75 nominal, a relationship in which fi is less than 1 . 

6. A method as claimed in Claim 4, in which 
the modification effected in step Ig) of the factor 
K is carried out by taking the factor K as K 
nominal 



80 



85 



90 



7. A method as claimed in Claim 4. in which 
the modification effected in step (g) of the factor 
K is carried out by taking the factor K as K 
nominal. 

I-. 

8. A method as claimed in Claim 4. in which 
the modification effected in step Ig) of the factor 
K is carried out by taking the factor K as K 
nomine.4 



in which 



n+f(L,-L,)). 



f{L,.L,) 



is a function of the difference between Lj and L,. 
9. A device for cjriYtng out the method as 
95 claimed in any one of Claims 1 to 8, comprising a 
gauge for measuring the level of the fluid 
contained in the lank, means for measuring the 
amount of fluid withdrawn from the tank per unit 
of lime, and means for determining a first value 

100 L,lnl related to the amount of fluid present in the 
tank 10 depending on the information (L^) 
representative of the fluid level, weighting this 
information on the basis a first value L,(n- 1 ) 
previously determined, further for determining a 

1 05 second value related to the amount of fluid 

present in the tank from the measured flow rate 
to which is applied said factor K. further for 
comparing the two values L,. L,. further for 
resetting the second value to equal L. when 

no the difference between the values L,. L, exceeds 
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a given threshold, and for modifying the factor K 
depending on the drfference existing between the 
first and second vaiues and Lj. 

1 0. A device as claimed in Claim 9. in which 
5 said means for measuring the amount of fluid 

withdrawn from the tank is a flowmeter. 

1 1. A device as claimed in Claim 9 or 10, 
comprising means for storing the first and second 
values L,, and means for displaying the second 

1 0 value to represent the amount of fuel present In 
the tank. 

1 2. A device as claimed In any one of Claims 9 
to 1 1Jn which said means for determining said 
first value Is arranged to calculated the value 



15 



(p-1) Lo 
L,(n)« L,(n-1)+ 

P P 



in which: 

L,(n) represents the first value L, at the 
moment (n(, 

L,(n~l ) represents the first value L, at the 
20 moment (n— 1 ), 

1^ represents a datum proportional to the 
information representative of the level of 
fluid delivered by the gauge (11). and 

p represents any whole number. 

25 1 3. A method as claimed In claim 1 and 

substantially as lierein described, for measuring 
the amount of fluid contained in a tank. 

14. A device for measuring the amount of fluid 
contained in a tank, which device is substantially 

30 as herein described with reference to the 
accompanying drawings. 
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METHOD AND APPARATUS FOR MONITORING 
OPERATIONAL PERFORMANCE OP FLUID STORAGE SYSTEMS 
Background of the Inventiion 
5 Th6 Invention relates to monitoring the 

Operational performance of fluid storage systems* 

Large quantities of liquids and similar materials 
are often stored in bulk storage containers or tanks, 
which may be located above-ground, partially above- 
10 ground, or completely below ground. Such containers or 
tanks are generally connected by piping to flow-meters or 
dispensers* 

For example, underground storage tanks (UST's) 
and, occasionally, above--ground storage tanks (AST's) are 

15 used to store petroleum products and fuel to be dispensed 
at automobile service stations, trucking terminals, 
automobile rental outlets, and similar operations through 
gasoline, diesel, or kerosene dispensing pimps. Fuel 
product is generally delivered to such facilities by a 

20 gravity drop from a compartment in a wheeled transport 
means such as a fuel delivery truck. AST's or UST's are 
often located at central distribution locations so that 
product can be subsequently withdrawn from the tank 
system to be transported for delivery to a variety of 

25 such facilities. A distribution location with UST's or 
AST's may receive deliveries of product from, e.g., a 
pipeline spur, wheeled transport, a barge, or a rail car. 

Direct observation of the operating condition of 
such tanks and storage containers is difficult or 

30 impossible. The vaurious methods for identifying the 

amount of product in tank systems have varying levels of 
accuracy, repeatability, and performance. Moreover, the 
accuracy of devices which measure the amount of product 
dispensed from the storage containers and tanks differs 

35 greatly, and may or may not be temperature compensated. 
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The aaount of product actually delivered to the tank 
system is often measured inaccurately and, frequently, 
not at all. Rather, the owner of operator of the tank or 
vessel usually records the invoiced amount of product 
5 delivered as the actual amount introduced to the tank 
system, without having any meems of confirming whether 
the invoiced amount of product delivered is correct. 

Consequently, effective management of such 
facilities is complicated by the numerous errors in the 
10 various measuring devices and procedures used to 
establish a baseline for management, planning and 
decisionmaking. Effective management requires the 
following: 

1. Accurate measurement of the volume stored in the 
15 system. 

2. Accurate determination of the volume dispensed 
from the system. 

3. Accurate determination of the amount of product 
introduced into the system. 

20 4- Identification of volumes added to or removed from 

the tank system which are not otherwise recorded. 

5. Rapid identification of leakage from the tank 
system • 

6. Continuous monitoring and diagnosis of the 

25 operating performance of all of the component 

measuring devices of the system. 

7. Continuous analysis of sales data to predict 
demands of product from the system. 

a. Determination of optimal reorder times and 

30 quantities as a function of ordering, 

transportation, holding, and penalty costs in 
order to minimize total costs of operation and/ or 
to maximize profits. 
' Trsiditionally, these functions were performed 
35 prudely or, in many cases, not at all. Volume 
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measurements were, and in many instances still are, based 
on inqperfect knowledge of the geometry, dimensions, and 
conf igiiration of the storage vessel. Also, dispensing 
meters are frequently aiscalibrated* This is true even 
5 when tank systems are regulated, due to the breadth of 
tolerance permitted for individual sales as related to 
total tank volume. For example, deliveries from the 
delivery vehicle are almost always unmetered, additions 
of product from defueling vehicles are typically 

10 undocumented, and theft of the product is not uncommon. 

Leakage of product has, in recent years, assumed a 
dimension far in excess of the mere loss of the product. 
E:nvironmental damage can, and frequently does, expose the 
operator to very large liabilities from third party 

15 litigation in addition to U.S. Environmental Protection 
Agency (EPA) -mandated remediation which can cost in the 
range of hundreds of thousands of dollars. The EPA's 
requirements for leeOc detection are set forth in EPA Pub. 
No. 510-K-95-003, Str^^jght T^lk On TftnKss Leg^K Detectjop 

20 Methods For Petro leum Underground Storage Tanks and 
Piping (July 1991) , which is incorporated herein by 
reference. 

To address these concerns, Statistical Inventory 
Reconciliation (SIR) was developed. The SIR method 

25 consists of a computer-based procedure which identifies 
all of the sources of error noted above by statistical 
analysis of the various and unique patterns that are 
introduced into the inventory data and, in particular, 
into the cumulative veuriances in the data when viewed as 

30 functions of product height, sales volumes, and time. 
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Pymn?irY invention 
The present invention relates to an autoBatic SIR 
system that may continuously and automatically collect 
data from completely above-ground, partially above- 
5 ground, and completely below ground containers for 

statistical analysis. The invention addresses a variety 
of physical, business, operational and environmental 
issues associated with the bulk storage of liquids or 
pour able solids. 

10 The present invention is an application of SIR 

that greatly enhances the ability to manage a facility 
effectively. It provides the means to characterize 
exactly the geometry, dimensions, and configuration of 
the storage vessel, identify overages and shortages in 

15 deliveries and imexplained additions and removals of 
product, and provide an accurate assessment of overall 
dispensing meter calibration. In addition, by accounting 
for such discrepancies, the present invention permits 
identification of leakage at rates less than .1 gallon 

20 per hour in all of its estimates to any prescribed 
tolerance. By increasing the number of measurements 
taken, the estimates can be derived at any desired level 
of tolctrance. 

The method of the present invention makes no 

25 assumptions as to the precision of any of the measuring 
devices used in various system configurations. Precision 
and calibration accuracies are derived from the data 
alone. Also, it is not assumed that the tank system is 
leak free; the leak status of the system is determined 

30 from the data alone. 

The method derives tank geometry, dimensions, and 
configuration, and their impact on the totality of 
cumulative inventory variances, as a function of product 
height in the tank. Correctness of dispensing meter 

35 calibration is verified in a similar maimer by testing 



wo 97/46855 



PCT/US97/09505 



- 5 - 

for randomness of cumulative variances as a function of 
varying sales volumes. Having confirmed that such 
remaining residual variances are random, reflecting only 
the inherent random noise of the measurement devices, the 
5 present method analyzes departures of the cumulative 
variance from the bounds determined by the calculated 
random noise level. All calculations as to the volumes 
added, removed, metered or leaking are based upon 
extended successive, simultaneous observations of meter 

10 and gauge readings. The number of observations 

incorporated in each such calculation is determined by 
computing confidence bands for the parameters of interest 
and extending data collection as necessary to achieve 
predetermined tolerances. 

15 For example, the method of the present invention 

is capable of distinguishing between continuous losses 
consistent with leakage and one-time unexplained removals 
of the fluid product from the tank. The method may be 
used to ensure the accuracy of computed delivery voliames, 

20 which are determined and reported with confidence 

boundaries calculated for estimated delivered quantities. 

The method can also be used to control and monitor 
the accuracy of purchase costs of fluids such as 
petroleum which are delivered to tanks. For exeunple, 

25 motor fuel retailers may be charged by wholesalers for 
either net or gross volumes purported to have been 
delivered. A determination that purchase charges are 
appropriate thus rec|uires frequent simultaneous readings 
of sales, tank volumes and temperatures, which can be 

30 accomplished using the method of the present invention. 

To accomplish these goals, the present invention 
involves estimating changes in product volume in a tank 
based on multiple data points and their respective likely 
errors* measured continuously over a period of time. A 

35 software program is used to inclement an algorithm that 
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employs concepts from matrix theory and mathematical 
statistics. The algorithm includes generating the 
product of a matrix and its transpose by successive 
additions of partial products of partitions of the matrix 
5 and their corresponding transposed matrix partitions to 
minimize the storage requirements of the data collected. 
The compressed matrix data constitutes a complete and 
sufficient statistic for the parameters of interest. The 
algorithm thus permits the accumulation and storage of a 

10 large amount of data in a condensed form without 

sacrificing statistically useful information, to obtain a 
statistically significant result with the required 
accuracy and reliability. 

Thus, one object of the present invention is to 

15 determine the accuracy and consistency of devices used to 
measure volume of product added to, removed from, and 
present in a fluid storage system. 

Another object of the invention is to identify and 
quantify additions of material to the system, but not 

20 recorded as such, and volumes of product removed from the 
system which are not registered by measuring devices or 
otherwise recorded. 

Another object of the invention is to discriminate 
between discrete one-time unrecorded removals of product 

25 from the system and continuous losses consistent with 
leakage. 

Another object of the invention is to identify and 
provide early warning of product leakage from all parts 
of the system, extending from the fill point to the point 
30 of discharge, and to confirm the validity of resulting 
leakage warnings. 

Another object of the invention is to determine 
secular, seasonal trends and repetitive special demands 
to provide short and long term estimates for demand of 
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the product, and to provide optimal reorder quantities 
and a delivery schedule for the system. 

A further object of this invention is to 
accomplish all of the foregoing in a fully automated 
5 system that recjuires no human intervention, other them as 
an option available to the operator to enter quantities 
of material reportedly delivered for comparison with 
those computed. 

In general, in one aspect, the invention features 

10 a method of monitoring a fluid storage and dispensing 
system* The system has measurement apparatus for 
measuring a volume associated with the system • A 
plurality of measurement data is collected from the 
measurement apparatus in a form readable by a computer 

15 and stored in a compressed matrix format in a computer 
memory. The compressed matrix foirmat is statistically 
analyzed to determine operational monitoring information. 

Implementation of the invention also may include 
one or more of the following features. The statistically 

20 analyzing step includes calculating error data resulting 
from the measurement apparatus. The method may also 
include determining the presence of operational defects 
in the system, monitoring the accuracy of the measurement 
apparatus, determining the volume of fluid in the system, 

25 and determining whether a cpiantity of fluid removed from 
the system is caused by a leak in the system. The method 
also may include delivering a warning that a leak has 
been detected in the system. 

The collecting step may be performed while the 

30 system is operating, may be performed continuously at 
periodic intervals, and may be performed automatically 
under the control of the computer. The method may 
include querying the measurement apparatus under the 
control of the computer. 
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The method may include measuring the temperature 
of the fluid and calculating the volume based on the 
temperature of the volume in the system. The calculating 
step may include determining a correction value based on 
5 a weighted average of the temperature of the volxime 
measured at a plurality of locations within the system. 

The system may include a plurality of tanks, 
including an underground storage tank, an above-ground 
storage tank, or a partially above-ground storage tank* 

10 The storing step may include generating the 

compressed matrix format as a product of a data matrix 
and the transpose of the data matrix. Further, the 
product may be formed by addition of partial products of 
each of a plurality of partitions of the data matrix with 

IS the transpose of each of the partitions. 

The measurement apparatus may include a volumetric 
gauge, a dispensing apparatus and a sales recording 
device. The measurement data or the compressed matrix 
format may be transmitted to a host processor to perform 

20 the statistical analysis. The collecting step also may 
include estimating an initial volume of fluid in the 
system. 

In general, in another aspect, the invention 
features a method of monitoring a fluid storage and 

25 dispensing system. The system includes measurement 
apparatus for measuring a volume associated with the 
system. A plurality of measurement data is collected 
from the measurement apparatus in a form read2Q)le by a 
computer and stored in a compressed matrix format in a 

30 computer memory. The compressed matrix format is 

statistically analyzed, to determine whether a quemtity 
of fluid removed from the system is caused by a leak in 
the system. A warning that a leak has been detected in 
the system is delivered. 
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In general, in another aspect, the invention 
features an apparatus for monitoring a fluid storage and 
dispensing systeai. Measurement apparatus measure a 
volume associated with the system. A computer includes a 
5 processing means for collecting a plurality of 

measurement data from the measurement apparatus and a 
memory for storing the plurality of measurement data in a 
compressed matrix format. The processing means performs 
statistical analysis of the compressed matrix format to 
10 determine operational monitoring information. 

Implementation of the invention also may include 
one or more of the following features. The system may 
include a plurality of tanks, including an underground 
storage tank, an above-ground storage tank, or a 
15 partially above-ground storage tamk. 

The measiirement apparatus may include a volumetric 
gauge, a dispensing appcuratus and a sales recording 
device. The volumetric gauge may include a probe 
disposed in the system, and the probe may be a 
20 magnetostrictive tank probe. The dispensing apparatus 
may include a totalizer. The sales recording device may 
simulate operation of a point of sale terminal. Also, 
the measxirement data may be collected simultaneously from 
the various measurement apparatus. 

The system may include a temperature sensor 
disposed in the system in contact with the volume for 
obtaining a temperature measurement in the volume. 
Alternately, the system may include a plurality of 
temperature sensors disposed at different locations 
within the system, each of the sensors being in contact 
with the volume for obtaining a plurality of temperattire 
measurements of the volume. The temperature measurements 
may be used to calculate the volume. 

The measurement data obtained from the measurement 
apparatus may be transmitted to the computer. A host 
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processor remote froa the computer may perform the 
statistical analysis. 

The compressed matrix format may be generated as a 
product of a data matrix and the transpose of the data 
5 matrix. Further, the product may be formed by addition 
of partial products of each of a plurality of partitions 
of the data matrix with the transpose of each of the 
partitions • 

In general, in another aspect, the invention 

10 features a method of monitoring a fluid storage and 
dispensing system. The system has a plurality of 
measurement apparatus for measuring a volume associated 
with the system* Measurement data is simultaneously 
collected from the plurality of measurement apparatus in 

15 a form readable by a computer to determine a change in 
the volume. The collecting step is repeated to obtain a 
plurality of measurement data from the plurality of 
measurement apparatus* The plurality of measurement data 
is stored in a compressed matrix format in a computer 

20 memory, and the compressed matrix format is statistically 
analyzed to determine operational monitoring information* 

Implementation of the invention may also include 
the following feature. The method may include estimating 
an initial value of the volume during the analyzing step. 

25 In general, in another aspect, the invention 

features a method of monitoring a fluid storage and 
dispensing system. The system has measurement apparatus 
for measuring a volume associated with the system and a 
plurality of temperaitiire sensing devices located at 

30 different heights in the system. The volume has a height 
in the system. A plurality of volume measurement data is 
collected from the measurement apparatus in a form 
readable by a computer. The volume measurement data is 
^adjusted based on temperature measurements taken from 

35 those of the plurality of temperature sensing devices at 
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a height below the height of the volume in the system. 
The plurality of volume measurement data is stored in a 
compressed matrix format in a computer memory, and the 
compressed matrix format is statistically analyzed to 
5 determine operational monitoring information. 

In general, in another aspect, the invention 
features a method of determining a volume associated with 
a fluid storage and dispensing system. The volume has a 
height in the system, and the system has measurement 

10 apparatus for measuring the height. A plurality of 

height measurement data is collected from the measurement 
apparatus in a form readable by a computer. The 
plurality of height measurement data is stored in a 
compressed matrix format in a computer memory. 

15 Regression analysis is performed using the compressed 

matrix format to calculate the volume associated with the 
system. 

Implementation of the invention may also include 
one or more of the following features. The collecting 
20 step may be performed each time a portion of the volume 
is dispensed from the system. The collecting step may 
not be performed when the fluid is being added to the 
system. 

In general, in another aspect, the invention 
25 features an apparatus for determining a volume associated 
with a fluid storage and dispensing system. The volume 
has a height in the system, and measurement apparatus 
measiures the height of the volume. A computer has a 
processing means for collecting a plurality of height 
30 measurement data from the measurement apparatus and a 
memory for storing the pliurality of height measurement 
data in a compressed matrix format. The processing means 
performs regression analysis of the compressed matrix 
^format' to determine the volume associated with the 
35 system. 
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In general, in another aspect, the invention 
features a method of determining a plurality of voltunes, 
each of the voltmes associated with one of a plurality of 
fluid storage and dispensing systems. Each of the 
5 volumes has a height in its associated system, and each 
of the systems has measiurement apparatus for measuring 
the height for each of the volximes. A plurality of 
height measurement data is collected from the measurement 
apparatus of each of the plurality of systems in a form 
10 readable by a computer. The plurality of height 

measurement data is stored in a compressed matrix format 
in a computer memory. Regression analysis is performed 
using the compressed matrix format to calculate the 
volumes associated with the systems. 
15 Other features and advantages of the invention 

will become apparent from the following detailed 
description, and from the claims. 

Brief Description of the Drawings 
Fig. 1 is a schematic diagram of a facility 
20 including an underground tank storage system. 

Figs. 2, 3 and 4 are a portion of the Mathcad 
computer code used to perform the data compression 
algorithm. 

Figs. 5, 6 and 7 are a block diagram of the steps 
2S performed during routine operation of the algorithm of 
the present invention. 

Fig. 8 is a block diagram of the steps performed 
during the data deletion operation of the algorithm of 
the present invention. 
30 Pig. 9 is a block diagram of the steps performed 

during the delivery calculation operation of the 
algorithm of the present invention. 

Fig. 10 is a schematic diagram of a data 
'acquisition and transmission network that may be used in 
35 conjunction with the present invention. 
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Fig. 11 is a schematic diagreua of a facility 
including an above-ground tank storage systea. 

Fig. 12 is a schematic diagram of a facility 
including a partially above-ground tank storage system. 
5 Description of the Preferred Embodim^niis 

The method and apparatus described herein applies 
to UST's, AST's or any type of storage tank. The product 
stored in the tank may be any fluid, including dry 
particles that flow in the manner of a fluid. 

10 Fig. 1 shows a UST facility lo, illustrated as an 

automobile service station. Facility 10 includes a 
series of UST's 12, 14, 16 which may store the same or 
different types of liquid fuel product 18. Volumetric 
tank gauges 20, 22, 24 in each tank measure the height of 

15 product 18 in the tank. Submersible pumps 26, 28, 30 in 
each tank pump product 18 to one of dispensing pumps 32, 
34 through piping lines 36, 38, 40. Alternately, 
facility 10 may be an AST facility with above-ground tank 
1000, as shown in Fig. 11, or a facility with a partially 

20 above-ground tank 1010, as shown in Fig. 12. 

Tank gauges 20, 22, 24 are moxinted in tsmks 12, 
14, 16. The tank gauges may consist of or be based on 
magnetostrictive tank probes or other sensing 
technologies. In the case of magnetostrictive 

25 technology, two floats 42, 44 surround each probe, e.g., 
gauge 20 in tank 12. One float 42 floats on the upper 
siirface of product 18 in tank 12, and the other float 44 
floats on the interface of product 18 with any water or 
other foreign material collected at the bottom of tank 

30 12. Tank gauge 20 determines the distance between floats 
42, 44 to obtain the height of product 18 in tank 12. 
Tank gauge 20 also contains temperature sensors 46, 48, 
50 spaced along its length to monitor the. temperature of 
product 18 at various depth levels. 
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Each of the dispensing punps 32, 34 consists of a 
totalizer or flov meter 52, 54 disposed in a housing 56, 
58 to measure the volume of product 18 dispensed through 
hoses 60, 62 and nozzles 64, 66. To operate dispensing 
5 pump 32, nozzle 64 is removed from housing 56, which 

actuates dispensing pump 32 and causes product 18 to flov 
through hose 60 due to the pumping action of submersible 
pumps 26, 28, 30. A value stored in totalizer 52 is 
incremented as fuel is dispensed through hose 60. Upon 

10 completion of the transaction, nozzle 64 is replaced in 
housing 56, thereby turning off dispensing pump 32 and 
discontinuing the action of submersible pumps 26, 28, 30 
and totalizer 52. 

Transactions are recorded electronically by 

IS software in a sales recording device 71 connected to 

totalizers 52, 54 of dispensing pumps 32, 34. Totalizers 
52, 54 in dispensing pumps 32, 34 are connected to sales 
recording device 71 by means of communications and power 
supply vires 78, 80. 

20 Sales recording device 71 contains softvare 

capable of emulating the functions of a point of sale 
(PCS) terminal associated vith fuel sales made at 
facility 10. PCS emulation softvare in sales recording 
device 71 functions on the basis of read only commands to 

25 eliminate the possibility of conflict with control 
commands from a POS terminal employed by facility 10. 
Alternative data acquisition systems can result in 
destruction of credit card sales records, inadvertently 
shutting down the entire system, and/or causing 

30 electrical interference in the pump linXs. 

Tank gauges 20, 22, 24 are connected to a tank 
monitor 82 by means of communications and power supply 
wires 84, 86, 88 or communicate data through radio 
* frequency transmission. Tank monitor 82 converts raw 
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data obtained from tamk gauges 20, 22, 24 into a fora 
usable by a computer. 

A computer 70 contains a processor 72 capable of 
running various computer softweure applications and a 
5 memory 74. Tank monitor 82 and sales recording device 71 
are electrically connected to computer 70 to relay 
totalizer values, product height and temperature data to 
computer 70. Software executable by processor 72 of 
computer 70 is capaUdle of querying tank monitor 82 and 

10 sales recording device 71 to obtain measurement data at 
selected time intervals. The data is continuously 
evaluated as it is collected and is stored in memory 74 
of computer 70 for later retrieval and detailed analysis. 
Alternatively, computer 70 may communicate with a host 

15 processor 90 at a remote location. The continuous 

evaluations or detailed analysis may then be performed by 
host processor 90, which may be faster or more efficient 
than computer 70. 

As an example, computer 70 may be a personal 

20 computer or any other proprietary microprocessor-based 
unit. Computer 70 may capture data automatically through 
direct-*connect serial interfaces with tank monitor 82 and 
sales recording device 71, or by manual operator keypad 
entry. Computer 70 communicates with equipment at 

25 facility 10 through four programmable serial 

communication ports, such as RS*232 communication ports. 

Computer 70 may, e.g., store tank dimensions and 
product characteristics, and concurrent time and date 
data along with the measurement data. Computer 70 may be 

30 used to produce error and analysis reports as calculated 
by the software. It may also have alarm event-- initiated 
capabilities, such as when a leak is detected in any of 
the tanks. Such a computer system can accommodate 
^facility and customer specific requirements while 
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maintaining conplete compatibility with other system 
components • 

The SIR method Involves reconciling volume data 
obtained from tank monitor 82 and volume data obtained 
5 from sales records. Sales transactions may be detected 
in a number of ways, including an electronic signal 
emitted from totalizers 52, 54, by voltage sensing of 
control relays on pump dispensers 22, 24, or by 
observation of product removal using tank gauges 20, 22, 
10 24. 

It is essential that the measurements used to 
obtain these two types of data are made simultaneously. 
The SIR method of the present invention collects and 
analyzes observations of sales volumes and tank volumes 

15 which are derived simultaneously. Failure to collect 
both types of data simultaneously would bias estimates 
derived from separate volume measurements. 

The SIR method properly accounts for the effects 
of temperature, pressure and specific gravity. In 

20 addition, product from two or more tanks may be blended, 
such as to achieve varying petroleiim octane levels at 
pump dispensers 32, 34. When different fluid products 
are blended, the tamks are treated as one unit, and an 
additional parameter is introduced to determine the 

25 actual blend percentages. 

Data concerning the physical characteristics of 
the tank configurations and the accuracy of the various 
gauges and metering devices is collected during 
installation and a set-up phase of operation of facility 

30 10 to create a basis for subsequent statistical analysis. 
Information is then continuously collected so that the 
statistical analysis of SIR can be performed by computer 
70 or host processor 90. 

Several procedures are used either singly or in 

35 combination to obtain the volume observations. First, 
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where the system configuration provides for deteraining 
whether hoses and dispensers associated with a given tank 
are active, the system is queried on a minute-by-minute 
basis, or on the basis of another predetermined time 
5 interval, to determine the status of the dispensers. 
When all of the dispensers are idle, the values from 
totalizers 52, 54, the tank volumes (i.e. product heights 
in the tanks) and temperatures are recorded. 

Second, submersible pumps 26, 28, 30 are checked 
10 to determine on/off status. When it is determined that 
the pumps are turned off, the values from totalizers 52, 
54 are read, and tank volumes and temperatures are 
recorded • 

Third, softwaure algorithms used by computer 70 

15 detect and measure leads and/or lags between the 

recording of sales events and corresponding gauge and 
meter readings. When leads or lags are encotintered and 
constitute a physical characteristic of the data 
measurement and recording system, constrained 

20 optimization, rather than unconstrained optimization, may 
be used to determine pareuneter estimates. Lagrange 
multipliers are one example of such a constrained 
optimization method. 

The method of the present invention is capable of 

25 providing dyneunic monitoring of system performance. For 
example, the leak detection function is carried out 
continuously while normal operations, e.g., removals and 
deliveries, are taking place. To detect leaks 
dyneunically, the software is programmed to detect when 

30 sales or delivery events occur and to calculate the 
volumes of product removed or added as a result of such 
activities. Thus, dynamic testing does not require that 
the system be dormant and addresses the entire system 
from the point of filling to the point of dispensing. 
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The SIR method of the present invention also 
distinguishes between one-tine removals and continuous 
losses consistent with lealcage. The integrity or lea3c- 
free status of the system is not assumed a priori* 
5 Instead, the individual and unic[ue characteristic pattern 
induced by each form of error when viewed along the 
separate dimensions of time, product height and sales 
volume are used to identify and quantify the errors. The 
method may also be tised to detect and quantify 
10 undocumented removals, e.g., theft or additions of 
product . 

Further, the overall system is self diagnosing in 
that it determines from the data the maximum degrees of 
reliability and precision of which a particular operating 

15 configuration is capable at any given time, as well as 
the degree of calibration accxiracy. 

In particular, product height in the tanks and 
temperature are measured continuously at, e.g., one- 
minute intervals. Height and gross volumes are converted 

20 to net volumes at, e.g., 60**? or 15*C, using the 
algorithms described below. Sales recorded by the 
totalizers 52, 54 are extracted and stored in memory 74 
at times coincident with readings from tank gauges 20, 
22, 24. If the dispensing system is capable of 

25 transmitting a signal indicating whether or not any or 
all individual hoses are active, that information is also 
stored in memory 74 coincident with taking gauge and 
meter readings. 

The method of the present invention is designed to 

30 achieve the maximum accuracy possible within the 
limitations imposed by the inherent random and 
irreducible noise in the various measuring devices 
incorporated. It utilizes multiple measurements over 
extended time periods to identify and quantify systematic 

35 and repeatable effects in the instrumentation and thereby 
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correct for such effects using the known physical 
characteristics of the devices. The system makes no a 
priori assumptions as to the accuracy of the devices used 
to ncasure product volume in the tank, to measure volumes 
5 removed, or as to the accviracy of volumes reported to 
have been delivered into the system. 

The resulting volumetric calculations are 
independent of the physical characteristics of the tank 
configuration and the various measuring devices which may 
be incorporated in the system. The results do not rely 
on input entered externally by the operator or from 
diagnostics internal to the measuring devices used. 
Instead, the output produced by the softwcure which 
analyzes the measured data depends only the patterns 
induced in inventory data produced by the tank gauges and 
measuring devices and, in particular, the cumulative 
variances that result when the various input values are 
combined. 

Various error patterns which the measuring devices 
can induce and the effects of temperature, tank geometry, 
and orientation on cumulative variances are derived from 
empirical analysis of real*world inventory data. The 
system^ s software synthesizes the output measurements of 
the various devices based on known characteristics 
derived from the empirical data. Thus, the software is 
capable of identifying measurement errors caused by the 
measuring devices and simultaneously compensating for the 
effects of those errors. 

Gauges can be systematically inaccurate in two 
ways. The height of the product in the tank can be 
incorrect, and the height to volume conversion algorithms 
may not reflect accurately the true dimensions of the 
tank or its orientation in the ground. The latter may be 
the result of incorrect measurements or an inappropriate 
conversion algorithm* 
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The presence of such systeiaatic effects and their 
nature may be established by examining the pattern of 
inventory variances as a function of product height. 
£rrors of this kind induce patterns which repeat 
5 themselves as the tank is filled and emptied. If the 
tank length is incorrect, a linear pattern is induced. 
If product height is in error, a curvilinear pattern 
results reflecting the varying volumes in different cross 
sections of a cylindrical tank. Tilt along the length of 

10 the tank induces a sinusoidal pattern symmetrical about 
the mid<*height of the tank. Absent such errors, the 
pattern will be purely random, reflecting only the 
inherent noise of the measxxring devices. The absence of 
randomness and the presence of a systematic pattern 

15 serves to identify the presence of systematic error. The 
pattern of a departure from random and its extent 
determines the source and extent of the effects and the 
means necessary to correct them. 

Dispensing errors, unlike volume measuring errors, 

20 are independent of product height, but are sensitive to 
the volume of product dispensed. The nature and extent 
of dispensing errors can be established by examining 
inventory variances as a function of sales volume. As in 
the case of volume measurements, in the absence of 

25 systematic errors, variances as a function of sales 

volume will be random. The form and extent of departures 
from randomness serve to determine the source and extent 
of the errors and provide for their removal. 

Leakage from the system creates a continuous 

30 downward trend in the cumulative variance when viewed as 
a function of time. By contrast, one-time additions and 
removals of product cause significant upward or downward 
translations of the cumulative variance which remain 
^permanently in the record and do not introduce a 

35 continuous trend. Leakage is distinguishable from tank 
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gauging errors when vieved as a function of product 
height because the pattern does not repeat as the tank is 
filled and emptied. If product is leaking from the 
system, a series of pairallel translations in the 
5 ciuQulative variance is generated, each shifted by the 
volume of product lost between deliveries. 

The accuracy of measurements taken from the 
various components of the system determines the accuracy 
achievable in any one individual observation. Since the 

10 lezOc rate is computed from a series of successive 

observations, however, the minimum detectable leak rate 
can be reduced to any desired magnitude by increasing the 
number of successive observations recorded. Thus, the 
system can serve as a final verification for leakage 

IS indications obtained by other methods. 

At the conclusion of an initial set up period of 
data collection including one or more delivery and sales 
cycles, the collected measurement data is analyzed by 
regression analysis. The initial set-up regression is 

20 used to derive tank dimensions and orientation, 

individual meter calibrations and secular trends. A 
confidence level value p is computed at the .01 level of 
significance to determine the minimum leak rate 
detectable by the system, and the residual variance is 

25 computed to provide the current noise level of the 
system. 

The regression is performed according to the 
following equation: 

stAR,L,T) ^a-il S a^a„i £l>^ - Eti Ls f^Bi (1) 

where : 

30 st^(R,L,T) = Volume in gallons derived from the ith 

gauge reading in inches in a cylindrical 
tank with or without hemispherical end 
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caps with radius R, length L, and tilt: 
over its length of T inches. 

a « Initial inventory in gallons, ii^ich is 
to be estimated. 

5 Sa^^j B Sales volume recorded on the Irth 

totalizer. 

ay^ = Fraction of sales volume recorded on the 
irth totalizer actually removed from the 
tank, which is to be estimated. 

10 Dj = Volume of the jth delivery. 

Et£ = Elapsed time since initiation until the 
ith gauge reading is recorded. 

Ls s Constant gain or loss in product per 
unit of time. 

15 Bj = Volume of product added (e.g. delivery) 

or removed during some discrete time 
interval prior to or during observation 
period j. 



0 if j < i 



All of the parameters are estimated simultaneously 
20 using least square estimation procedures. The R and T 
parameters are derived numerically, but the other 
parameters are derived analytically. 

Further, all of the parcuneter s , including the 
initial inventory, are estimated simultaneously. The 
25 initial volume must be estimated from all succeeding 

data, even if the tank is initially empty, otherwise the 
initial gauge reading and its conversion to gallons is 
assigned a credibility not assumed for all succeeding 
readings. Also, in a great majority of applications, the 
30 initial inventory in an already existing and operating 
-system. is not accurately known. 

Initial inventory estimation is vital in 
determining the geometry of the tank. When tank 
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geometry, tank orientation, or temk product height 
measxirement depart from the values obtained from nominal 
sources, all gauge and meter measurements are affected. 
It is practically impossible to detect the errors induced 
5 in the gauge measurements and correct for them unless the 
estimation of the initial inventory is made coincident 
with the estimation of the values of the other 
parameters. 

The estimate of the parameters are based on the 
10 totality of the data collected. This means, e.g., that 
the estimate of leak rate Ls is determined from a linear 
trend including all of the data collected, not merely at 
one end of the reconciliation period. Likewise, 
estimates of tank dimensions and orientation are derived 
15 from their overall contribution to reduction in residual 
variance, as opposed to a sale by sale analysis of tank 
segments. 

The volume st£(R,L,T) is derived from the product 
height measurement by multiplying the constant area of 
20 tank segments of height h (in inches) by tank length L. 
The volxime in gallons of product in a horizontal 
cylindrical tank of radius R is given by: 

Vol = ^ [i?2cos-^ ' IR ' h) {2Rh - h^) ^] 



In the case of a tilted tank, the area of the 
segments vary with position along the length of the 
25 tilted tank, and the volume is determined by integrating 
over the length L. Such integration does not result in a 
^closed form because the cross sections are not circular, 
and a numerical integration would severely limit the 
frequency of observations. Instead, in this application 



I 
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the l:anlc is treated as lying horizontally and the product 
is considered tilted, to derive an equivalent volume* 
This integration yields the closed form: 

Vol = JL^ [ (z - Dsin-^ (2z . 2 - 4 (2z - z^) M * w 
231 3 — 

The integrand is evaluated between the normalized product 
5 heights in inches, hu/R and hl/R, at the lower and higher 
ends of the tilted tank, respectively. It is standard 
industry practice to Install tanks on an incline to 
divert water and sludge away from the submersible pumps. 
Tank tilt is identified from the pattern it 

10 induces in the record of cumulative variances as a 
fxinction of product height. It is compensated for by 
fitting the correct mathematical form for height to 
volume conversions in a tilted temk to the cumulative 
variance calculated by the method of least squares. This 

15 is done simultaneously with estimation of the initial 
inventory . 

Tank length L and radius R are established by 
equating the first partial derivatives of the sum of 
squared cumulative variance with respect to length and 

20 radius and determining the values which minimize the sum 
of squared variances. Simultaneous estimation of initial 
inventory is also required when estimating tank length L 
and radius R. 

Errors in measurement of the product height h in 

25 the tank are characterized by curvilinear patterns 

induced by height to volume conversions in the cumulative 
variance for a cylindrical container when heights are 
transposed upward or downward. Such errors also are 
compensated for by minimizing the sum of squared 

30 ^cumulative variances with respect to increments or 

decrements to measured product height. This estimation 
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also requires simultaneous estimation of the initial 
inventory of the tank. 

In general, the accuracy of the estimates of the 
tank: dimensions, tank orientation and height measurements 
5 is confirmed by observing that the cumulative variemces 
of each derived value as a function of nominal product 
height are random and display no systematic influence or 
effects. 

Dispenser totalizer calibration is continuously 
10 monitored and evaluated by minimizing the sum of squared 
cumulative variances with respect to multiplicative 
constants associated with individual reported cumulative 
sales volumes from all pump dispensers associated with a 
particular tank system* This eliminates the need for 
manual verification of meter calibration. 

In particular, gauge performance is continuously 
monitored to identify gauge malfunctions or degradation 
in gauge performance. Monitoring of gauge performance is 
independent of diagnostics which are internal to the 
measuring device. Diagnoses of problems are based only 
on their impact on the cumulative inventory variances 
which are continuously monitored by the software. 

If the gauge fails to record changes in product 
height when the dispensers register sales, an increase in 
cumulative variances approximately equal to sales volume 
is observed; this effect can be identified by the 
monitoring software and a warning of gauge malfunction 
generated to the operator. 

However, observation of the gauge registering 
product height change, but with a time lag after sales 
are recorded, may be a feature of normal gauge 
performance. Such normal gauge performance is identified 
by repeated positive increments in cumulative variances 
as sales are completed with subsequent return of the 
ctunulative variance to normal bounds. When such gauge 
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function is determined to be the normal operating 
characteristic of a particular system, constrained 
optimization with lagged variables is introduced into the 
software, otherwise, the gauge's performance is reported 
5 as a malfunction* 

Finally, temperatures in the tank are monitored to 
detect changes that are excessive for the time intervals 
between observations. Erratic temperatiire readings are 
deleted, and may indicate gauge malfunction. 

10 The software computes actual, rather than nominal, 

delivered quantities and requires no input by the system 
operator. The operator may choose to input into the 
system the nominal delivery quantity indicated by the 
delivery invoice, along with the temperature and 

15 coefficient of expansion of the product at the point of 
pick-up. The software will then compute overages or 
shortages between the nominal and actual quantities 
delivered, as well as the overages or shortages caused by 
temperature-* induced variations in the transport of the 

20 product to the facility and in the subsequent mixing of 
the delivered product with that resident in the tank. 

Delivery is identified by the software when a 
positive cumulative variance is observed which exceeds 
the system noise level and is not succeeded by a return 

25 to normal variance boxinds. Delivered quantities are 
computed by estimating the volume increases they induce 
in multiple, successive observations. The required 
number of successive observations is determined as that 
sufficient to generate a confidence width which is within 

30 a predetermined tolerance. The system of the present 
invention is capable of accounting for sales conducted 
during delivery and for noise introduced by post delivery 
turbulence in the tank. 

One-time unaccounted for removals or additions to 

35 the tank are computed in the same manner. Deliveries are 
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distinguished from such events by computing the rate of 
input, which in the case of nonaal gravity delivery 
should exceed 100 gallons per minute. Other modes of 
delivery, e.g. pipeline delivery into above ground tanks, 
5 are identified by incorporating their known delivery 
rates. 

Leakage from the system is identified by a 
continuous linear negative trend in the data which 
exceeds the computed minimum detectable leak rate after 

10 all of the various error phenomena described above have 
been identified and compensated for. This calculation 
deals with the totality of the data obtained by 
constantly monitoring kno%im removals and is not 
restricted to observations made only when the system is 

15 dormant. It is also independent of any single data 

reconciliation calculation in that trends throughout all 
of the data are evaluated. 

All calculations concerning volumes are made on 
the basis of net volumes, according to the following 

20 definitions: 

Net Volume 
in Tank 

where: 

t 

25 

CE 

and 

Net Sales 
30 Volume 

where t^^ and tj are temperatures measured by the tank 
"gauge at the beginning and ending of a sale transaction, 
respectively. Deliveries are computed in net gallons. 



Gauge Volume (1 - (t-60)CE} 



Measured temperature in degrees 
Fahrenheit (if centigrade, the term in 
parentheses becomes (t-15) ) . 

Coefficient of expansion. 



Metered Sale (1 - ( t^ ^ t^ - 60) CE) 

2 
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but are converted to gross quantities if required, based 
on external information input by the system operator, as 
follows: 

GT « Gross gallons on invoice at the 
5 originating terminal. 

NT « Net gallons on invoice at the terminal. 
tT = Temperature at the terminal. 
CE ■= Coefficient of expansion. 

The program also records: 

10 tA = Ambient temperature in the tank prior to 

delivery. 

tF = Temperature in the tank at the 
conclusion of delivery. 

The following value is computed: 

15 tS ^ Temperature of the product in the 

delivery vehicle at the facility at the 
beginning of delivery. 



tF + Mh (tF - tA) 
NVD 



20 where: 



NVD - Actual net volume delivered, previously 
computed. 

NVA s Net volume in the storage tank at the 
start of delivery. 

25 NS » Net overage(+) (underage (-)) in 

delivery. 

NT-NVD 

GVD =» Gross volume delivered. 

NVD (1 + (tF - 60) CE) 

30 GVS = Gross volume in the transport vehicle at 

the facility prior to delivery. 

NVD (1 + (tS - 60) CE) 
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GSH = Shrinkage due to mixing in the tank. 
GVS - GVD 

GVT « Actual gross volume in the transport 
vehicle at the facility. 

5 » NVD (1 + (tT - 60) CE) 

GST = Shrinkage during transit to the 
facility. 

GVT - GVS 

GOS « Gross overage (+) (underage (-) ) adjusted 
10 for temperature effects 

GT - GVD + GST + GSM 

Calculations of voliunes actually delivered are 
based on multiple observations of the balance of measured 
tank volumes and cumulative sales. This method recfuires 
15 frequent simultaneous observations of sales and in-tank 
volumes (i.e. product heights) and temperatures. 

The volume of product in a tank is derived by 
measuring the height of the product and using the 
geometry of the tank, which is assumed to be known, to 
20 compute the corresponding volume. In many instances, 
tank dimensions vary substantially from assumed design 
dimensions. Regulatory specifications permit up to 10* 
variation in length and diameter of cylindrical tanks. 

Tank orientation can also cause complications in 
25 the calculations. The volume corresponding to a measured 
height varies substantially when the tank is tilted away 
from horizontal or rolled away from vertical. 

Further, temks may also fail to conform to a known 
geometry either through faulty manufacture or 
30 installation, or may suffer significant deformation 
during the course of operations. For example, many 
-fiberglass tanks sag or bend along their length. 

In addition, installed temks are typically 
inaccessible, and difficult to measure. Thus, it is 
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necessary to confirm the accuracy of height to volume 
conversions from generated inventory data and to identify 
and correct discrepancies where they exist. 

The foregoing problems are compounded when two or 
5 more tanks are manifolded together* Manifolded tanks are 
joined together by piping systems and serve common 
dispensers. Thus, sales quantities from manifolded tanks 
constitute withdrawals from all tanks in the manifolded 
system, but not necessarily in equal quantities. Product 

10 heights typically vary from tank to tank, but tank 

geometries, dimensions and orientation may also vary so 
that a procedure for correcting height to volume 
conversion errors for a single tank will not apply. 

The different factors which influence inventory 

15 data manifest themselves in distinct ways which 

facilitate their identification and correction. These 
factors are most easily identified by examination of 
their effects on cumulative departures of actual measured 
inventory from a theoretical or book value when viewed 

20 across a variety of dimensions. In particular, one-time 
undocumented physical additions or removals of product, 
e.g. over or under deliveries and pilferage, are 
evidenced by an addition or subtraction of a constant 
quantity from the cumulative variance at the time of 

25 occurrence and all subsequent observations. Continuous 
loss of product accumulating over time, e.g. leakage, is 
evidenced by a loss trend over time. Continuous loss of 
product varying proportionally with sales value, such a 
line leeOc or meter miscalibration, may be determined by 

30 identifying a constant negative trend that is c\mulative 
only over periods where delivery lines are pressurized. 

A pattern of gains or losses, or both, recurring 
cyclically as the tamk is successively filled and emptied 
"with no long term gain or loss of product, is the pattern 

35 associated with height to volume conversion error. The 
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pattern is cyclical because the error source is identical 
in each cycle as the tank is filled or emptied. It is 
distinguishable from the other patterns in that it 
retraces the same path without the translation which 
5 would occur if physical loss or gain of product were 
taking place. 

This problem is most readily diagnosed by 
analyzing cumulative variance as a function of product 
height. If the variances are random with no evidence of 
systematic effects, height to volume conversions may be 
assumed to be correct. If not, the form of the induced 
pattern indicates the nature of the conversion error. 
Thus, an error in tank length induces a linear pattern, 
an error in tank tilt induces a sinusoidal pattern, and a 
constant error in tank height measurement induces an arc- 
like pattern. When other sources of loss or gain are 
present, the conversion error patterns remain, but are 
translated in each succeeding filling/emptying cycle to 
reflect the physical loss of product which has occurred 
during that cycle. Thus, confusion between conversion 
errors and other effects can be eliminated. 

Sales readings and product height measurements 
must be made simultaneously. Since the number of 
observations in any one sales cycle is typically too few 
to generate a conversion table of sufficient detail to be 
of practical use, subsequent sales cycles and their 
corresponding deliveries must be incorporated. If, 
however, deliveries are umnetered and are used to 
approximate the volume (as is the standard industry 
practice), significant inconsistencies are introduced. 
If an overage or shortage occurs dxuring delivery, then 
all subsequent sales volumes correspond to tank cross 
sections which have been shifted upward or downward from 
their predecessors. Averaging or statistical treatment 
cannot overcome this deficiency since there is no means 
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of Icnovlng without metering whether « by how much, and in 
what direction the data has been shifted. 

The procedure of the present invention may include 
determining if height to volume conversion error is a 
5 problem. If the error is a problem, then the system must 
determine the nature of the problem, e.g. tank 
dimensions, tank orientation, height measurement or 
unknown tank geometry, and whether the conversion problem 
is compounded by other gains and losses. If leakage is 

10 suspected, an on-site leak detection investigation is 
undertaken. In no leakage is indicated, and one or all 
of tank dimensions, tank orientation and height 
measurement are problems, new conversion factors are 
calculated and confirmed using the diagnostic procedures 

15 described herein. 

If unknown tank geometry or manifolded systems are 
encountered, the exact current percentage of metered 
sales actually dispensed from each dispenser is 
determined by physical measurement. A high order 

20 polynomial using a variable of measured product height is 
used to convert height to volume. The parameters of the 
polynomial are derived from the differences between 
measured product height corresponding to the beginning 
and ending of sales events which do not overlap 

25 deliveries. 

For a single tank, actual dispensed quantities are 
regressed using a polynomial based on the differences in 
measured product height before and after individual 
sales, subject to the constraint that when the polynomial 

30 is evaluated at a height equal to tank diameter, the 
result is the total tank volume. Observations which 
include delivery events are discarded. 

ASal^i « ^li^i-i - ^i) * *2(**^f.i - h^i) ♦ • • • 

35 Vol « a^d + agd^ + . . . ♦ a^d" 
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A fifth order polynomial has proven adequate in 

most cases* Residual analysis may be used to determine 

adec[uacy of the polynoiaial in the presence of severe tank 

distortions, and higher order polynomials may be 

s introduced as necessary. The number of observations 

required is determined by estimating a confidence bound 

around the resulting polynomial with a width adequate for 

the desired resolution. Thus, 

ASale^^ - Actual dispensed volume in period 
10 i. 

h^ = Product height upon conclusion of 

ASalej^. 

= Product height prior to 

commencement of ASale^ and 
15 after completion of ASale^.^^, 

d ^ Diameter of tank. 

Vol " Total volume of tank. 

The converted volume for height h is then given by: 

Vol (h) « a^h + ajh^ + . . • + a„h" 

20 The omission of a constant term in the regression implies 
that 

Vol (h) = 0 when h « 0 

This ensures that the polynomial derived from the height 
differences is well defined. 

25 For manifolded systems, actual sales are regressed 

simultaneously on individual polynomials based on the 
various height differences in the several tanks which 
correspond to a particular sales volume, subject to the 
constraint that each polynomial evaluated at the 

30 corresponding tank diauneter yields the total volume of 
that tank. 
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where: 



5 ASale^ = Actual Sales voltune in period !• 

h^.^j = Height of product in tank j after 

completion of Asale^^^^ and prior to 
cononencing Asale^^. 

j = 1, 2, ... m 

10 h^j « Height of product in tank j after 

completion of ASale^. 

a a Niimber of tanks manifolded. 

Volume conversion for the tn measured heights, h^^, h2f • . 
. h„ in the total system is: 

15 m n 

Vol (h^, h2, . . . h^) = £ r a^i h^^ 

i-1 J-1 

where: 

- Height of product measured in the 
20 ith tank in the manifold. 

Delivery inaccuracies have no impact on this 
calculation since all observations made during deliveries 
are discarded. Height changes are related only to the 
corresponding volumes dispensed. 

25 Prior determination of actual quantities 

dispensed, as opposed to metered quantities, ensures that 
the only remaining source of error is random measurement 
error. Regression is designed to accommodate random 
error of this kind and to facilitate inferences when 

30 errors are present. 

With respect to temperature, the temperature of 
■product delivered into a tank system almost invzuriably 
differs from the temperature of the product already in 
the tank. Its addition has the effect of escpanding or 
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contracting the voltme of the combined product. This 
change in volume can create the appearance of incorrect 
dimensions of the height to volume conversion, appear as 
leeJcage where none exists, or it can mask the existence 
5 of actual leakage* 

It is therefore preferable, and frequently 
essential, that all volumes, sales, deliveries and 
product in storage be converted to a common temperattire 
prior to analysis. Typically 60*P (IS^'C) is chosen as 
10 the standard. The conversion is accomplished as follovs: 
Net Volume = Gross Volume (1 - (t«-60)CE) 



where: 



Measured product temperature in degrees 
Fahrenheit. 



15 CE = Coefficient of expansion. 

As above, all calculations are in net gallons of product. 

A complication to the calculation may occur if the 
tank gauges 20, 22, 24 used to measure product volume are 
designed for static or dormant mode tank testing. Such 

20 tank gauges detect leakage when the tank is taken out of 
service. In this case, product volume changes due to 
temperature changes during the course of a test must be 
accounted for. 

Further, as shown in Fig. 1, temperature sensors 

25 46, 48, 50 are located at different heights in tank 12. 
If the level of product falls below a given temperature 
sensor, the corresponding weighted temperature 
measurement is dropped from the average temperature 
calculation, and a temperature jump and corresponding 

30 volume change may be observed when the net volume is 
calculated using the new weighted average of 
temperatures. If uncorrected, such repeated jumps in the 
data would preclude further analysis of the data for leak 
detection or the generation of height to volume 

35 conversions. 
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The system of the present invention may be used to 
overcome these temperature related problems. Using the 
following definition, 



NDB^i » The net cumulative variance in the 
inventory data at observation 



then. 



NDB„ = a(l - {to - 60) CE) 
N 

-r Sa^Cl - (t| - 60) CE) 
10 i-1 

-VK(l-(tN - 60) CE) 

Where: 

a = Gross initial inventory. 

to = Temperature of initial product volume. 

15 t^ = Temperature of product at observation i. 

Sa^^ - Gross volume sold in period i. 

V|| e Measured gross volume in tank at period 
N. 

CE = Coefficient of expansion. 

20 Absent random error or leakage, and assuming no 
deliveries of product, then 



NDB|| » 0 



and 



wo 97/46855 



FCT/US97/09505 



- 37 - 

a(l - (to • 60) CE) 



N 

S Sa^ (1 - (t^ - 60) CE) 



5 « (1 - (t^ - 60) CE) 

Therefore, if a temperature jximp to temperature t* occurs 
at an observation N+1, then 

NDBr+I = Vn(1 - (t„ - 60) CE) 

- Sa^+i (1 - (t* - 60) CE) 
" - Vh^I (1 - (t* - 60) CE) 

V^d - (tn - 60) CE) 

" Sa^+i (1 - (t* - 60) CE) 

- (V„ - Saj,+i) (I - (t* - 60) CE) 
V„(t*-t^)CE 

^5 When this final quantity NDB^^.^ is added to the 

volume where the transition occurs between temperature 
sensors, and all subsequent volumes, the effect of the 
transition is eliminated, and analysis proceeds as it 
would where individual temperature readings are 

20 available. 

A large number of variables must be estimated by 
the software to implement the SIR system of the present 
invention. For example, as many as forty hoses and 
independent totalizers per tank system, as well as 

25 deliveries numbering four or more per day must be 

accommodated. Thus, a very large volume of data must be 
accumulated, encompassing a substantial spread of sales 
volumes from each totalizer for both the set-up analysis 
■and subsequent routine monitoring. To accommodate this 

30 volume of data within current or conceivable future 
practical computer memory capabilities, the algorithm 
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implemented by the software utilizes a matrix formulation 
which invoices the property of a sufficient statistic to 
reduce the memory requirement. 

The calculations used to determine the various 
5 error, loss trend and delivery estimates have the form: 

B = (x^x)*^x^y 

MSB » fy-xP)^(y-yP) 
m+1 

= (x^x)'^ MSB 



xo where : 



B « Column vector of m parameters to be 
estimated • 



X = Matrix of pareuneter coefficients. 

y B Column vector of independent variables. 

15 MSB — Mean squared error. 

» Variance covariance matrix of parameter 
estimates . 

The values contained in vector y comprise tank 
gauge readings. The entries in matrix x are measured 

20 sales volumes, time, and other constant values. The 

parameters of vector B which are to be evaluated include 
the initial volume of the system and subsequent volume 
changes, including delivery amounts. 

For example, if observations are recorded every 

25 minute, as many as 1440 rows in the x matrix and the y 
vector may be recorded. It would clearly be impractical 
to accumulate and store data in that form over an 
extended period of time. Instead, data compression 
techniques are applied so that only a manageable amount 

30 of data need be stored. 

The algorithm utilizes the property that if an n x 
m matrix A is partitioned into two submatrices, B and C, 
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where B is an i x m matrix and C is a j x a matrix, such 
that i + j « n, then 

C^C « A^A + B^B 
For example, 
5 Ti 2 

if C = / I 2 then A 



15 



10 



A^A = 



B^B = 



Thus, A^A + B^B 




[1 2 1 and B =[2 ll 

■ [::] 



[2 .2 1 2] 



1 2 

1 2 

2 1 
1 2 



[::,] 



At the conclusion of each 24 hour or other period, 
20 only x^x and x^ computed and stored. The matrix x 

has the form of a square n x n matrix. Further, the 
aggregates of observations for different periods are 
additive, since two square matrices having n x n 
dimensions may be added. Thus the total data storage 
25 requirement for each period is determined only by the 
square of the number of parameters of interest. 

The system is able to accommodate virtually 
unlimited numbers of observations by this method of data 
compression. Without this capability, the system would 
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not have the storage capacity to accurately and 
siaultaneously estimate the numbers of parameters which 
are required to perform a statistically significant 
calculation. This data compression method also allows 
5 for processing the data at the facility or for 

transmitting the data to a host computer for periodic 
analysis. Pigs. 2, 3 and 4 show the Mathcad computer 
code used to perform the data compression algorithm. 
Furthermore, (x^x)*^x^y is a complete and 

10 sufficient statistic for B. No statistically useful 
information is lost in the compression. The overall 
procedure is, therefore, unlimited by memory. The only 
limitation remaining is the precision available in the 
computer system used. 

15 The software performs SIR analysis, including 

inventory estimation and leak detection, using the above 
equation in the following form: 

y = (Stk^ange) " (^^^range) 

20 Where : 

range - 1 . . • (number of observations) 

meters * 1 . • • (number of dispensers) 



and 



^range ^ Column of I'S. 



25 "^range * Cumulative time in minutes. 



# Q \ <meterB> 

\^ranga/ 



Cumulative sales for an individual 
dispenser in gallons. 



^range ~ Cumulative deliveries. 
30 S'^^anga ^ Tank stick reading in gallons. 

To estimate the initial inventory, the matrix x 
includes a column of unitary values. To estimate loss 
-trends I the matrix x includes a column containing 
cximulative times of measurement and cumulative sales. 
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The values of B, MSB and are then calculated, 
producing the following result for the vector B: 
Bi = Estimated initial inventory. 

B2 - Loss trend. 
^ ^2-«'meter8 " Individual meter error. 

B is the vector containing the parameter 
estimates, namely beginning inventory, meter calibrations 
and loss rate. The loss rate estimate is in the second 
row (n»2) • 

10 s^ is the variance covariance matrix of the parameter 

estimates. Thus, = (5^22)^'^ is the standard deviation 
of the loss rate estimate. Finally, the minimal 
detectable leak is defined as t^Sjs. where t^ is the (l-o) 
percentile of the Student's t distribution. 

15 The software performs delivery calculations using 

the equation in the following form: 

y = (Stk^ange) 

20 Where: 

range = 1 . • . (number of records) 

and 

^rang© = Column of I's. 
^rango Cumulative time in minutes. 

25 ^rango * Cumulative sales in gallons. 

^ranga ^ where T^^nga is ^®ss than delivery time 

and 

1 where Ty^j,— is greater than or equal 

to 

30 delivery time. 

S^^range " T^hk Stick reading in gallons. 

The values of B, MSB and are then calculated, 
producing the following result for the vector B: 
= Estimated initial inventory. 
35 B2 = Loss trend. 
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B3 • Meter error. 

B4 » Estimated delivery amount. 

is the variance covariance matrix of the 
estimates. Thus^ S44 = (S^^^)^^^ is the standard deviation 
5 of B4, the delivery volume estimate. The delivery 

tolerance is B4±t^S44, where t^ is the (l-a> percentile of 
the Student's t distribution. Delivery tolerances can be 
reduced to any desired value by increasing the number of 
observations used in the calculation. 

10 The SIR analysis used by the method of the present 

invention involves computing and comparing cumulative 
variances. When the initial set-up is complete, computed 
trend and meter calibrations are used to project forward 
an expected cumulative variance, that is, the expected 

15 value of the difference between gauge readings and 

computed inventory. Actual cumulative variances are then 
computed from all subsequent gauge and meter readings and 
compared to the expected variance. 

Figs. 5, 6 and 7 show the routine operation 

20 procedure 100 followed by the software to perform this 
analysis. Data from the set-up of the system and the 
most recent analysis is entered into the program at step 
102. The data entered includes the tank type, tank 
dimensions, tank tilt, meter calibrations, mean square 

25 error atnd calculated trends. At step 104, three 

variables established as counters, Counterl, Counter2, 
and Counter3, are set at zero. The measurement data from 
the system itself is entered at step 106, namely the 
readings from the dispenser totalizers, the product 

30 height and the product temperature. 

The software computes the gross vol\ime of the 
product, the most recent gross volume and the sales as 
measured by the individual dispensers at step 108. The 
software further manipulates the data at step 110 by 

35 converting all gross volumes to net volumes, computing 
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observation to observation variance, and cosaputing 
cumulative variance. The sign of the cumulative variance 
is recorded at step 112. 

The program proceeds on the basis of the 
5 cumulative variance and the value of Counterl in steps 
114, 120, 124, 128, 132, 136 and 140. Depending on the 
cumulative variance and the value of Counterl, the 
program analyzes the collected data at step 118 if it is 
a final observation (step 116), deletes. the collected 

10 data (steps 122 and 134), performs the analysis for a 

delivery (step 126) (see below) , or reads new data (steps 
116, 122, 130, 134, 138 and 142) upon updating the value 
of Counterl and other computational variables (i.e. 
index, sign index and sign) • In some cases, collected 

15 data is deleted (steps 122 and 134). 

Upon computing the rate loss at step 144, the 
program reads new data at step 146 if the rate loss is 
not greater than or equal to, e.g., .2 gallon per hour, 
otherwise it computes the trend of the data at step 148. 

20 If at step 150 it is determined that the trend is greater 
than .2 gallon per hour, a warning is issued at step 156. 
In either case, the software continues to read and 
analyze the data at steps 152, 154, 158 and 160 until the 
last observation. 

25 The operation of deleting data 170 is shown in 

detail in Fig. 8. After performing similar analyses at 
steps 172, 174, 176, 180 and 184, using the indices and 
the values of the calculated standard deviations as in 
the routine operation procedure described above, the 

30 values of the counters are updated and new data is read 
at steps 178, 182, 186 and 188. Data is deleted in 
accordance with steps 178, 186 and 188. 

Finally, Pig. 9 shows the delivery calculation 190 
in detiaiil. After determining that the cumulative 

35 variance is greater than a predetermined value (three 
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Standard deviations) at step 192, the program determines 
whether the variance Is greater than, e.g., lOO gallons 
per minute (step 194). If so, the delivery Is recorded 
and the amount of delivered is determined at steps 202 
5 through 220. 

If there is a delivery in progress (step 202) , 
data is read until a negative observation to obsexrvatlon 
variance is observed (step 204) . The variance is 
monitored until the turbulence in the tauik subsides (step 

10 206) • Thirty observations are read (step 208) , and all 
observations from 15 minutes before the delivery until 
the end of the turbulence observations are deleted (step 
210) . An indicator variable is introduced with the 
turbulence observation, froB which regression commences 

IS (step 212). The confidence bound on the indicator is 

computed (step 214). If the confidence bound is within a 
predetermined tolerance, the volume of the delivery is 
reported within the confidence bounds (step 220) ; 
otherwise, additional observations are added, and the 

20 confidence bound is recomputed (step 218). 

If the variance between data measurements is less 
than 100 gallons per minute, the software determines 
whether the gauges are inoperative and reports them as 
being inoperative (step 198), or proceeds as in the 

25 routine operation procedure according to step 200 (in 
which there is a negative variance) depending on whether 
the observation exceeds a predetermined value (within one 
standard deviation) at step 196. 

In general, if observed variances are within three 

30 standard deviations or other predetermined tolerance of 
the expected value, the data is stored for future 
analysis* When cumulative variance exceeds three 
standard deviations or other predetermined tolerance, 
different software programs are executed depending on the 

35 nature and magnitude of the departure. 
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If within ten (or other predetermined) successive 
observations after the initial departure, the cumulative 
variance returns to within the tolerance range, all data 
froa and including ths initial departure and prior to the 
5 initial observation are deleted. The time extent and 
number of observations involved is recorded and stored 
for, e.g., a daily gauge performance report. 

If all ten (or other predetermined) successive 
observations remain outside the tolerance bound and the 

10 cumulative variances are of the same sign, a new trend 
line is initiated at the point of initial departure. 
After ten (or other predetermined) additional 
observations, a third trend line is initiated. If the 
increment to the overall trend estimated from the most 

15 recent observations is not significant, the most recent 
data is consolidated with the previous data and the 
process is repeated until such time, if ever, that the 
cxirrent trend increment is significant. 

If the departure is positive, the system checks 

20 whether the product is being dispensed and whether the 
gauge height fails to decrease, reflecting removal from 
tank. If so, the tank gauge is reported to be 
inoperative . 

If the gauge height is increasing, monitoring is 
25 continued as above until the most recent trend line 
returns to its original slope. Minute to minute 
variances are monitored to detect turbulence until the 
gauge values again return to within tolerance. All 
observations which occurred in the fifteen minutes prior 
30 to first positive departure \mtil the end of post 

delivery turbulence are deleted. An indicator variable 
is introduced at the first observation after post 
delivery turbulence. The system collects thirty 
additional observations and performs the regression from 
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the beginning of the period to determine the volume 
delivered. The volume delivered is then reported. 

If the depetrture is negative, the system proceeds 
as with delivery. If successive slope increments fail to 
show a return to the original slope, indicating 
continuing loss of product for a predetezi&ined period, 
typically one hour, and slope exceeds .2 gallon per hour, 
the system reports a warning that there is a continuous 
loss of product. If the loss rate is less than .2 gallon 
per hour but greater than the minimum detecteJsle leak, 
the system continues to monitor and recalculate the 
parameters, to be included in a daily operational report. 
If the incremental trend line shows a retxirn to the 
original trend, the system proceeds as with delivery, 
introduces an indicator variable, deletes data as 
necessary, and performs the regression to determine the 
volume of product removed. The system reports a one-*time 
removal of product. 

Referring to Fig* 10, the invention incorporates a 
data acquisition and transmission network (DAT Network) 
300 to completely automate the process of obtaining, 
capturing, transferring and processing product inventory 
data for use in product management, delivery scheduling 
and environmental compliance practices. DAT network 3O0 
includes on-site processors 302, 304 at the facilities 
306, 308 where the tanks are located, a customer host 
processor 310 and a central host processor 312, DAT 
network 300 links multiple remote facilities 306, 308 to 
central host processor 312, which performs the SIR 
analysis. The link may be accomplished indirectly 
through customer host processor 310, which itself is 
connected to a plurality of remote facilities 306, 308. 
Each of these processor elements is con^osed of 
"independently operating software and hardware systems 
which form the basis of a wide area network linked by 
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modems which transmit information electronically via the 
telephone network 314 using standard dial-up voice grade 
telephone lines. Examples of OAT networks are the 
TeleSIRA and ECCOSIRA systems developed by Warren Rogers 
5 Associates, Inc., Middletown, Rhode Island. 

DAT network 300 provides a uniform method of 
integrated management for the widest possible variation 
of underground and above-ground fuel storage, movement 
and measurement systems. On-site processors 302, 304 are 

10 capable of obtaining information from any electronic or 
mechanical control system, enabling DAT network 300 to 
accommodate facility configurations that are unique to 
each facility while presenting the information captured 
at remote facilities 306, 308 to customer host processor 

15 310 or central host processor 312 in a uniform format. 
On-line prbcessors 302, 304 obtain and capture 
product inventoiry data through the use of proprietary 
interfaces with external systems in use at remote 
facility 306, 308, such as tank gauges and sales 

20 recording devices. On-line processors 302, 304 transfer 
captured information daily, weekly or monthly through the 
public switched telephone network 314 to customer host 
processor 310 or central host processor 3X2 for use in 
inventory management, delivery scheduling and/or 

25 environmental compliance. On-site processors 3 02, 304 
may be, e.g., a touch-tone telephones acting as a sending 
units and Windows-based multi-line, voice prompt/ response 
PC's as the receiving units. On-site processors 302, 304 
may be designed to meet the specific needs of facilities 

30 306, 308 without requiring remote hardware at the 
facility in addition to that already present. 

In particular, each of on-site processors 302, 304 
typically is equipped with an alphanumeric keypad, a 
"character display, a power supply, four programmable 

35 serial communication ports, an internal auto-dial /auto- 
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answer (AD/AA) modem and a local printer port (for 
connection to a printer) . The keypad and display allow 
for operator configuration and manual entry of sales, 
delivery and tank level data* Use cf an AD/AA 2400 baud 
5 modem allows multiple on-site processor 302, 304 to share 
an existing voice grade telephone line by establishing 
communication windows to minimize attempted simultaneous 
use. Each of the programmable serial communication ports 
is independent, fully progrsunmable and governed by 

10 options selected at the facility or off -site through 
modem access. Finally, on-site processors 302, 304 can 
prompt the facility operator to enter missing or suspect 
entries when results are outside the expected range. 

The use of customer host processor 310, which is 

15 capeUt>le of receiving, storing and processing information 
from multiple on-site processors 302, 304, enables the 
management of a remote tank population from a single 
point of contact. A database of information created by 
customer host processor 310 is the basis for all higher 

20 level product management functions performed by DAT 
network 300. The database is also the basis for the 
environmental compliance analysis performed by central 
host processor 312. 

The use of central host processor 312, which is 

25 capable of receiving, storing and processing the 

information in the database created by customer host 
processor 310 for product management enables DAT network 
300 to achieve maximum results by utilizing the database 
for environmental compliance without additional remote 

30 facility information or communication. Central host 
processor 312 Is capable of transmitting a resulting 
database of the environmental analysis back to customer 
host processor 310 for printing and other customer 
record-keeping requirements. 
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The processor elements of DAT network 300 nay 
exhibit other useful operational characteristics. To 
prevent unauthorized access to DAT network 300, a 
security access code for dial-up data transfer functions 
5 is required. Under secured access, the baud rate, 
parity, stop bit parameters and communication protocol 
are determined at any of on-site processors 302, 304, 
customer host processor 310 or central host processor 
312. 

10 Another function of DAT network 300 is to monitor 

tank contents generally. DAT network 300 can be 
programmed to activate, e»g., an audible and visual alarm 
if the water level in the tank is too high (e.g., greater 
than 2 inches) , if the product level in the tank is too 

15 high (e.g., more than 90% of tank capacity) or too low 
(e.g., less than 10% of tank capacity, more product must 
be reordered, less than two days supply), and if a theft 
occurs (product level changes during quiet periods) . 
The system may be used to obtain valucO^le 

20 information other than inventory regulation and leak, 
detection. For example, the system may incorporate time 
series analysis routines, including Box Jenkins, moving 
average and exponential smoothing, to derive estimates of 
demand for the product which also incorporate temporal 

25 and seasonal trends and special events. 

The demand analysis may also be combined with 
additional inputs of holding costs, reorder costs, 
transportation costs and penalty costs for running out of 
stock. The system can include optimal inventory 

30 algorithms to determine optimal order quantities, reorder 
points and optimal delivery truck routing. Further, the 
system may incorporate multiechelon, optimal inventory 
procedures to accommodate combined wholesale and retail 
operations, such as with calculus-based optimization and 

35 linear, nonlinear and dynamic programming. 
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Other embodiiaents are within the scope of the 

claims. 
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What is claimed is: 

1. A method of monitoring a fluid storage and 



dispensing system, said system comprising measurement 
apparatus for measuring a volume associated with said 
5 system, said method comprising: 

collecting a plurality of measurement data from 
said measurement apparatus in a form readable by a 
computer; 

storing said plurality of measurement data in a 
10 compressed matrix format in a computer memory; and 

statistically analyzing said compressed matrix 
format to determine operational monitoring information. 

2. The method of claim 1 wherein said 
statistically analyzing step includes calculating error 

15 data resulting froia said measurement apparatus* 

3. The method of claim 1 further comprising 
determining the presence of operational defects in 

said system. 



20 



4. The method of claim 1 further comprising 
monitoring the accuracy of said measurement 



apparatus • 



5. The method of claim 1 further comprising 
determining the volume of fluid in said system. 



25 



6. The method of claim 1 fxirther comprising 
determining Aether a quantity of fluid removed 
from said system is caused by a leak in said system. 



7. The method of claim 6 further comprising 
delivering a warning that a leak has been detected 
in said system. 
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8. The method of claim 1 wherein said collecting 
step is performed while said system is operating. 

9. The method of claim 8 wherein said collecting 
step is performed continuously at periodic intervals. 

s 10. The method of claim 9 wherein said collecting 

step is performed automatically under the control of said 
computer. 

11. The method of claim 1 further comprising 
querying said measurement apparatus under the 
10 control of said computer. 



12. The method of claim 1 further comprising 
measuring the temperature of the fluid; and 
calculating said volume based on the temperature 
of said volume in said system. 



15 13. The method of claim 12 wherein said 

calculating includes determining a correction value based 
on a weighted average of the temperature of said volume 
measured at a plurality of locations within said system. 

14. The method of claim 1 wherein said system 
20 comprises a plurality of tanks. 

15. The method of claim 1 wherein said system 
comprises an underground storage tank. 

16. The method of claim 1 wherein said system 
comprises an above*ground storage tank. 

25 ' 17. The method of claim 1 wherein said system 

comprises a partially above*ground storage tank. 
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ls* The method of claim 1 wherein said storing 
step comprises generating said compressed matrix format 
as a product of a data matrix and the transpose of said 
data matrix. 

5 19. The method of claim 18 wherein said product 

is formed by addition of partial products of each of a 
plurality of partitions of said data matrix with the 
transpose of each said partition. 

20. The method of claim 1 wherein said 

10 measurement apparatus includes a volumetric gauge, a 
dispensing apparatus and a sales recording device. 

21. The method of claim 1 further comprising 
transmitting said measurement data to a host 

processor to perform said statistically analyzing step. 

15 22. The method of claim 1 further comprising 

transmitting said compressed matrix format to a 
host computer to perform said statistically analyzing 
step. 

23. The method of claim 1 wherein said collecting 
20 step comprises estimating an initial volume of fluid in 
said system. 



wo 97/46855 



PCTAJS97y09505 



- 54 - 

24. A method of monitoring a fluid storage and 
dispensing system, said system comprising measurement 
apparatus for measuring a volume associated with said 
system, said method comprising: 

5 collecting a plurality of measurement data from 

said measurement apparatus in a form readable by a 
computer; 

storing said plurality of measurement data in a 
compressed matrix format in a computer memory; 
10 statistically analyzing said compressed matrix 

format; 

determining whether a quantity of fluid removed 
from said system is caused by a leak in said system; and 
delivering a warning that a leak has been detected 
15 in said system. 

25. An apparatus for monitoring a fluid storage 
and dispensing system, said apparatus comprising: 

measurement apparatus for measuring a volume 
associated with said system; and 
20 a computer comprising a processing means for 

collecting a plurality of measurement data from said 
measurement apparatus and a memory for storing said 
plurality of measurement data in a compressed matrix 
format; 

25 wherein said processing means performs statistical 

analysis of said compressed matrix format to determine 
operational monitoring information. 

26. The apparatus of claim 25 wherein said system 
comprises a plurality of tanks. 

30 27. The apparatus of claim 25 wherein said system 

"comprises an underground storage tank. 
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28.. The apparatus of claim 25 wherein said system 
comprises an above-ground storage tank. 

29. The apparatus of claim 25 wherein said system 
comprises a partially above-ground storage tank. 

5 30. The apparatus of claim 25 wherein said 

measurement apparatus includes a voliimetric gauge, a 
dispensing apparatus and a sales recording device. 

31. The apparatus of claim 30 wherein said 
volumetric gauge comprises a probe disposed in said 

10 system. 

32. The apparatus of claim 31 wherein said probe 
is a magnetostrictive tank probe. 

33. The apparatus of claim 30 wherein said 
dispensing apparatus comprises a totalizer. 

15 34. The apparatus of claim 30 wherein said sales 

recording device simulates operation of a point of sale 
terminal . 

35. The apparatus of claim 25 further comprising 
a temperature sensor disposed within said system 

20 in contact with said volume for obtaining a temperature 
measurement of said volume. 

36. The apparatus of claim 35 wherein said 
temperature measurement is used to calculate said volume. 
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27. The apparatus of claim 25 further conprising 
a plurality of temperature sensors disposed at 
different locations within said system, each of said 
sensors being in contact with said volusa for obtaining a 
5 plurality of temperature measurements of said volume. 

38. The apparatus of claim 37 wherein said 
temperature measxirements are used to calculate said 
volume. 

39. The apparatus of claim 25 wherein said 

10 measurement data is collected simultaneously from said 
various measurement apparatus. 

40. The apparatus of claim 25 wherein said 
measurement data obtained from said measurement apparatus 
are transmitted to said computer. 

15 41. The apparatus of claim 25 further comprising 

a host processor remote from said computer to 
perform said statistical analysis. 

42. The apparatus of claim 25 wherein said 
compressed matrix format is generated as a product of a 

20 data matrix and the transpose of said data matrix. 

43. The apparatus of claim 42 wherein said 
product is formed by addition of partial products of each 
of a plurality of partitions of said data matrix with the 
transpose of each said partition. 
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44. A method of monitoring a fluid storage and 
dispensing system, said system comprising a plurality of 
measurement apparatus for measuring a volume associated 
with said system, said method comprising: 

5 simultaneously collecting measurement data from 

said plurality of measurement apparatus in a form 
readable by a computer to determine a change in said 
volume; 

repeating said collecting step to obtain a 
10 plurality of said measurement data from said plurality of 
measurement apparatus; 

storing said plurality of measurement data in a 
compressed matrix format in a computer memory; and 

statistically analyzing said compressed matrix 
15 format to determine operational monitoring information* 

45. The method of claim 44 further comprising 
estimating an initial value of said volume during 

said analyzing step. 
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46* A method of monitoring a fluid storage and 
dispensing system, said system comprising measurement 
apparatus for measuring a volume associated with said 
system and a plurality of temperature sensing devices 
s located at different heights in said system, said volume 
having a height in said system, said method comprising: 

collecting a plurality of volume measurement data 
from said measurement apparatus in a form readable by a 
computer; 

10 adjusting said volume measurement data based on 

temperature measurements taken from those of said 
plurality of temperature sensing devices at a height 
below the height of the volume in said system; 

storii^g said plurality of volume measurement data 
15 in a compressed matrix format in a computer memory; and 

statistically analyzing said compressed matrix 
format to determine operational monitoring information. 

47. A method of determining a volume associated 
with a fluid storage and dispensing system, said volume 

20 having a height in said system, said system comprising 
measurement apparatus for measuring said height, said 
method comprising: 

collecting a plurality of height measurement data 
from said measurement apparatus in a form readable by a 

25 computer; 

storing said plurality of height measurement data 
in a compressed matrix format in a computer memory; and 

performing regression analysis using said 
compressed matrix format to calculate said volume 
30 associated with said system. 

48. The method of claim 47 wherein said 
collecting step is performed each time a portion of said 
volume is dispensed from said system* 
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49. The method of claim 48 wherein said 
collecting step is not performed when said fluid is being 
added to said system. 

50. An apparatus for determining a volume 

s associated with a fluid storage and dispensing system, 
said volume having a height in said system, said 
apparatus comprising: 

measurement apparatus for measuring said height of 
said volume; and 
10 a computer comprising a processing means for 

collecting a plurality of height measurement data from 
said measurement apparatus and a memory for storing said 
plurality of height measurement data in a compressed 
matrix format; 

15 wherein said processing means performs regression 

analysis of said compressed matrix format to determine 
said volume associated with said system. 

51. A method of determining a plurality of 
volumes, each of said volumes associated with one of a 

20 plurality of fluid storage and dispensing systems, each 
of said volumes having a height in its associated system, 
each of said systems comprising measurement apparatus for 
measuring said height for each of said volumes, said 
method comprising: 

25 collecting a plurality of height measurement data 

from said measurement apparatus of each of said plurality 
of systems in a form readable by a computer; 

storing said plurality of height measurement data 
in a compressed matrix format in a computer memory; and 

30 performing regression analysis using said 

compressed matrix format to calculate said volumes 
associated with said systems. 
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